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ABSTRACT 


Recent power systems necessitate for maintaining a safe voltage stability as 
the number of problems such as contingencies and reactive power 
insufficiency are increasing. In this paper, installation and sizing of Flexible 
Alternating Current Transmission System (FACTS) devices have been 
introduced for solving the voltage stability problems under contingencies. 
The FACTS device to be used is Thyristor Controlled Series Capacitor 
(TCSC). Besides improving the voltage magnitude at all buses to a desired 
level, installation of TCSC at proper locations can minimize total 
transmission losses of the system. To conduct the sizing task, the newly 
developed Flower Pollination Algorithm (FPA) has been implemented as the 
engine for optimization. Through experimentation, the results proved that the 
proposed placement and sizing technique has successfully mitigated the 


Optimization voltage stability problems. In addition, the computation time for FPA’s 
Thyristor controlled series convergence was tolerable with optimum results. 
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1, INTRODUCTION 

The growing demand for electricity has made the utilities around the world as critical entity in 
energy market. To provide satisfactory services for consumers, a sophisticated scheme for maintaining 
voltage stability at a proper level is required. Voltage stability problem, especially voltage collapse, becomes 
a major issue due to the blackouts experienced by many countries for the past few years [1]. The term voltage 
stability is always referred to the ability of power system to maintain acceptable voltages as the system is 
subjected to a disturbance that may cause insufficiency of reactive power. Contingencies and failures of 
equipment are the types of disturbances occurred in power system which may lead to voltage instability [2]. 
Due to such problem, the power system may undergo the most serious malfunction: voltage collapse. 

Theoretically, insufficient reactive power support is the major cause for voltage instability where 
the voltage drop is uncontrollable [3]. There are variety of methods to mitigate the voltage instability 
problems and one of them is through the use of Flexible Alternating Current Transmission System (FACTS) 
devices. The Unified Power Flow Controller (UPFC), one of the FACTS devices, is effective as it can 
maintain reactive power requirement of the system in the event of large disturbance and faults [4]. Besides 
that, Thyristor Controlled Series Capacitor (TCSC) can be suitable device for improving and increasing the 
network power transfer capability [5]. Thus, it is understood that the FACTS devices improve the stability by 
controlling the level of reactive power in the system, which means through the injection and absorption of 
reactive power [6]-[7]. 
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To have a good installation of FACTS devices, proper placement (location for installation) and 
sizing (rating determination) should be conducted systematically as they promise the optimum improvement 
of voltage profiles and losses minimization [8]—[9]. While placement can be done through analytical methods 
such as sensitivity analysis and bus ranking, sizing is frequently conducted through optimization. The so- 
called ‘trial and error’ approach is rarely used when determining the size of FACTS devices as it does not 
offer optimal solution. It was said that proper selection of placement location can lead to efficient line flow 
control and maintain bus voltages at acceptable levels, thus enhancing the voltage dependability edges [10]. 
Implementing optimization approach necessitates for meta-heuristic algorithms like Genetic Algorithm (GA), 
Particle Sawrm Optimization (PSO), Ant Colony Optmization (ACO), Bacterial Foraging Algorithm (BFA) 
and the latest one Flower Pollination Algorithm (FPA). 

This paper proposes a method for optimal voltage stability improvement through TCSC and Flower 
Pollination Algorithm (FPA) considering contingency situations such as line outages. The study includes 
various comparative studies on the methods of TCSC placement as well as its effect on voltage profiles and 
transmission losses under various contingencies. Performance of FPA in the case study shall be assessed 
afterwards in terms of computation time for convergence and solution optimality. Lastly, all significant 
findings and analysis will be concluded in the last section of this paper. 


2. RESEARCH METHOD 
In this section, the fundamental of TCSC operation and proposed FPA algorithm are 
briefly explained. 


2.1. TCSC in Power System 

TCSC is a FACTS device installed in series with the transmission line of power system. 
Technically, TCSC consists of a series capacitor bank shunted by Thyristor Controlled Reactor (TCR). Such 
a connection provides a smooth variable series capacitive reactance, i.e. the control action taken on the 
reactance will give a corresponding reactive power compensation (either injection or absorption). 
The advantages of using TCSC in power system are voltage stability improvement, effective power flow 
control and damping power oscillations. Figure 1 shows the schematic diagram for TCSC connected along a 
transmission line. The whole circuitry which comprises of capacitor, inductor and two thyristors represents 
the basic TCSC circuit as in Figure | [11]. Hence, the combination of both capacitive and inductive reactance 
represents the TCSC’s equivalent reactance, X7c¢sc. 





Figure 1. Schematic diagram for TCSC 


In this paper, the range of TCSC reactance to be used is between -0.7 p.u. and 0.2 p.u. Thus, 
the developed FPA will search the best TCSC sizes between the range while ensuring no constraint violation 
during optimization. The algebraic summation between the line impedance and TCSC reactance will result in 
power flow control along the line, hence producing the desired compensation of reactive power. This control 
action will give a corresponding improvement on the bus voltages as well as losses reduction. Prior to that, 
however, the placement of TCSC will be conducted through the Maximum Loadability Identification (MLI) 
technique in which the sensitive lines may be identified. 


2.2. Proposed Algorithm for FPA 

Yang in [12] for the first time proposed an algorithm that mimics the pollination of flowers, 
known as the Flower Pollination Algorithm (FPA). The biological process of pollination involves pollinators, 
sometimes called as pollen-carrying agent, for spreading the pollens. Insects can be the example of 
pollinators. Basically, the FPA consists of two significant processes namely local and global pollination. 
Global pollination refers to biotic and cross-pollination whereby the pollinators travel in a path that obeys 
Levy flight. This means that the process necessitates for pollinators travelling at a long distance. On the other 
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hand, local pollination refers to abiotic and self-pollination. This process occurs when the pollen of the same 
flower fertilizes to be another flower. Usually, no pollinators such as insects involved in this process. 
The proposed problem formulation for this study is presented as follows: 


f(&}) = max {VY} (1) 
Where: 

Xi = [Xtcsc Xtescr = Xtesc -Xtesc (2) 
Subject to: 

—0.7 p.u.S Xresc S 0.2 p.u. (3) 

dy Pg = Pp + Pross (4) 

FVSI max < 0.9 (5) 


Where, X;’ is i-th solution (i.e. the pollen) at ¢-th iteration, X7csc is the TCSC reactance, N is the number of 
TCSC to be installed, Pg is the generator power, Pp is the total demand, P),,, is the total losses, FVSTjqy 1S a 
stability index poposed by Musirin [13] for constraint violation check, f(X;) is the fitness to be optimized and 
V, is the voltage magnitude at k-th bus. Hence, the objective function for FPA is to optimize the voltage 
profile at all buses. The proposed FPA algorithm in the context of the case study is presented as follows. 


Step 1: initialization 

At first, initialization of parameters is done by specifiying the values for FPA’s parameters 
heuristically. Next, a group of initial solutions, known as population, is generated randomly while satisfying 
all the constraints as in (3), (4) and (5). 
Step 2: fitness evaluation 

Later, all the randomly generated solutions are evaluated through load flow analysis. The aim is to 
improve the voltage magnitude at all buses as indicated in (1). 
Step 3: generate new pollens 

This is the step where the algorithm will decide either to use global or local pollination process in 
generating the new group of pollens, i.e. X;*’. First, in global pollination, flower pollen gametes are carried 
by pollinators such as insect, and pollen can travel over a long distance because insects can often fly and 
move in a much longer range. Therefore, it can be represented mathematically as: 


Xf** = XP + yLA)@* - X7) (6) 


Where, g is the current best solution found among all solutions at the current iteration. Here, y is a scaling 
factor to control step size. The scaling factor, y is set to 0.01. In addition, L(A) is the parameter that 
corresponds to the strength of pollination, which essentially is also the step size (i.e. s) with A is set to 1 x 10° 
° Since insects may move over a long distance with various distance steps, we can use Levy flight to imitate 
this characteristic efficiently; that is, to draw L(A) > 0 from a Levy distribution. This is given as follows: 


ar(a)sin(“*) 4 
a a glta ) 





L(A) = S>> So > 0 (7) 
Where, /(A) is the standard gamma function and this distribution is valid for large step size, s > 0. On the 
other hand, the local pollination is represented by the following equation: 


Xj** = Xj + e(Xj — X7) (8) 


Where, X; is pollen from different flower of the same plant species (i.e. random solution) and ¢ is a random 
number drawn from a uniform distribution between [0, 1]. 
Step 4: updating the population 

After the new pollens have been generated, they will be evaluated through step 2. Next, both the 
current and new population will be combined, sorted and discarded to maintain the original size of 
population. 
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Step 5: termination 
The convergence of the algorithm will be achived if all solutions have approximately equal fitness. 
Then, the optimization is terminated. Otherwise, step 2 to 5 will be repeated. The overall algorithm for FPA 


is Shown in Figure 2. 
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Figure 2. FPA’s overall algorithm 


3. RESULTS AND ANALYSIS 

There are three case studies focused here: (1) comparison between single and multiple installation of 
TCSC, (2) comparison between different methods of TCSC placement and (3) FPA performance. In this 
paper, four types of contingencies (or line outage) were specified and their notations are used as Cl, C2, C3 
and C4 to indicate their increasing severity. The IEEE 30-bus reliability test system (RTS) was used as the 
test system throughout this paper. 


3.1. Comparison between single and multiple installation of TCSC under Contingencies 

In this section, the effects of single and multiple installation of TCSC in power system are assessed 
based on minimum voltage magnitude and total losses at post optimization. The results are tabulated in 
Table | and the graphs as in Figure 3 and 4 were plotted based on the table. 
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Table 1. Improvement of voltage and total losses — single & multiple TCSC installation 


Before installation Single installation Multiple installation 
System Outage Voltage Losses Installed Voltage Losses Installed Voltage Losses 
state number (p.u.) (MW) line (p.u.) (MW) line (p.u.) (MW) 
Base - 0.9514 17.734 8 0.960 11.178 0.990 14.747 
Cl 1 line ().8984 64.041 8 0.960 42.532 6,8,10,17 0.980 27.845 
C2 2 lines ().8668 66.196 8 0.950 52.535 6,8,10,17 0.960 60.107 
C3 3 lines 0.8571 82.286 8 0.980 59.905 6,8,10,17 1.000 25.973 
C4 4 lines 0.8102 87.347 8 0.960 36.972 6,8,10,17 0.970 34.254 
1.05 100 


_0.95 Eo 8 so —_- 
J ay ™~ 


-— TAA 


20 


0.8 





cl c2 cs c4 . —_— cl cz cs c4 
Contingencies Contingencies 
Figure 3. Minimum voltage (in p.u.) for single and Figure 4. Total losses (in MW) for single and 
multiple installation of TCSC multiple installation of TCSC 


Based on Figure 3 and 4, with no TCSC installed to the system, the voltage reduces while the total 
losses increase as the contingency severity increases, that is, when the number of outage lines increases. 
The improvement on voltage profiles by single and multiple installation of TCSC is reasonable as the 
magnitudes are between 0.95 p.u. and 1.05 p.u. This can be seen in Figure 3. Through all contingencies, 
single installation promises an improvement to both voltage and total losses, but the multiple installation 
offers better results. For example, in Figure 4 the losses reduced by single installation at C3 is 59.905 MW, 
which is equivalent to 27 percent of reduction. A further reduction in losses is obtained when multiple 
installation is implemented, which is 25.973 MW or equivalent to 68 percent of losses reduction. It is not to 
say that multiple installation is the most recommended scheme, but if the cost of installation is not the issue, 
then it could be a favourable method among utilities or service providers. 


3.2. Comparison between Random and MLI-based Placement of TCSC 

Next, this section analyses the effects of two different methods of TCSC placement. The first 
method is by installing the TCSC at random locations in power systems, hence it is called as random 
placement. The second one selects the suitable locations based on a rank produced from a method called 
Maximum Loadability Identification (MLI), thus it is known as MLI-based placement. Article [14] explains 
the concept of MLI in the field of voltage stability assessment. The results are tabulated in Table 2 and the 
graphs as in Figure 5 and 6 were plotted based on the table. There is only a small improvement resulted from 
random placement as can be seen in Figure 5. At all levels of contingencies, the voltage magnitude is below 
0.95 p.u. This shows that inappropriate selection of TCSC locations does not promise any improvement to 
the system. 


Table 2. Improvement of Voltage and Total Losses — random and MLI-based Placement 


Before installation Random placement MLI-based placement 

System Outage Voltage Losses Installed Voltage Losses Installed Voltage Losses 
state number (p.u.) (MW) line (p.u.) (MW) line (p.u.) (MW) 
Base - 0.9514 17.734 8 0.960 11.178 8 0.960 11.178 
Cl 1 line 0.8984 64.041 6 0.910 63.801 13 0.960 62.313 
C2 2 lines 0.8668 66.196 16 0.870 66.198 8 0.950 35.590 
C3 3 lines 0.8571 82.286 10 0.880 81.766 12 0.960 30.527 
C4 4 lines 0.8102 87.347 17 0.820 87.341 2 0.950 47.624 
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Figure 5. Minimum voltage (in p.u.) forrandom and Figure 6. Total losses (in MW) for random and MLI- 
MLI-based placement of TCSC based placement of TCSC 


The MLI-based placement, on the other hand, has successfully maintained the voltage magnitude at 
least 0.95 p.u. at all levels of contingencies. Obviously, this method is consistent in promising an 
improvement to the system compared to random placement. The same goes to the losses whereby the random 
placement has failed to yield any obvious reduction. In contrast, the MLI results in effective minimization of 
losses especially at contingency C2, C3 and C4. The percentages of reduction at the three contingency levels 
are 19 percent, 63 percent and 45 percent respectively. Thus, it can be said that installation of TCSC 
necessitates for a proper method of placement in producing a better improvement to the system. 


3.3. FPA performance 

This section presents the analysis of FPA performance with regards to the two previous case studies. 
The parameter of FPA concerned here is the switch probability P,, which is used to control the preference of 
the algorithm, either global or local pollination, in producing new solutions. Two values were chosen 
intuitively to see the effect of the parameter on FPA performance, which are 0.5 and 0.8. The results are 
tabulated in Table 3 and the graphs as in Figure 7, 8 and 9 were plotted based on the table. Firstly, there is no 
obvious difference in voltage magnitude resulted from P, = 0.5 and P, = 0.8 as can be seen in Figure 7. 


Table 3. Effect of FPA’s switch Probability on Optimization Performance 


Single installation (P, = 0.5) Single installation (P, = 0.8) 
System Outage Installed Voltage Losses Time (s) Installed Voltage Losses Time (s) 
state number line (p.u.) (MW) line (p.u.) (MW) 
Base - 8 0.950 17.707 8.614 8 0.950 17725 9.945 
Cl 1 line 8 0.950 27.076 19.982 8 0.960 50.339 68.369 
C2 2 lines 8 0.950 27.206 67.936 8 0.950 35.879 66.474 
C3 3 lines 8 0.990 20.216 9.531 8 0.960 41.011 11.633 
C4 4 lines 8 0.960 29.809 10.218 8 0.960 41.771 10.546 
1.2 100 
a . 
0.8 
60 
0.6 
40 
0.4 
0.2 
; 0 


Figure 7. Minimum voltage (in p.u.) forP, =0.5 and _- Figure 8. Total losses (in MW) for for P, = 0.5 and P, 
P, = 0.8 — single installation = 0.8 — single installation 
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Figure 9. FPA’s computation time (in seconds) for P, = 0.5 and P, = 0.8 


The voltage magnitude for both values of P, is above 0.95 p.u. at all contingencies. However, this is 
different for losses minimization whereby the FPA did better if P, = 0.5. From Figure 8, the FPA with P, = 
0.5 was able to reduce more losses than that of P, = 0.8, and this is true for all levels of contingencies. 
Despite the voltage and losses at post-optimization, the required computation is another issue. A good 
optimization shall produce optimal solution with minimum computation time. Based on Figure 9, it is not to 
say that smaller P, is better than the larger one since both values required very close computation time, 
except for contingency Cl. This signifies that the switch probability P,, perhaps, has nothing to do with the 
speed of algorithm. It does affect the optimality of solution, but not the computation time. 


4. CONCLUSION 

In this paper, the optimization algorithm known as FPA has been developed for solving power 
system problem. The FPA was used to perform optimal sizing of TCSC for losses minimization and voltage 
magnitude improvement. The analysis of the results justified the ability of FPA in producing satisfactory 
voltage stability improvement at various conditions of system, including contingencies. Therefore, it is no 
doubt to conclude that the aims of this paper are successfully achieved. For future recommendation, the FPA 
can be hybridized with other meta-heuristic techniques to increase its performance in terms of convergence 
speed and optimality of solution. 
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